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Precede de localisation spatio-temporelle, de radio-mobiles en milieu urbain 


(57) 


Le procede consiste : 





a calculer (1) la reponse impulsionnelle du canal de 
propagation entre le radio-mobile et la station emet- 
trice/receptrice, 

a estimer (2) a partir de la r§ponse impulsionnelle 
les couples (angles, retard) du trajet de propagation 
direct et des multitrajets resultats de reflexions sur 
des obstacles situes dans ie champ de propagation 
des ondes entre la station de base et le radio- 
mobile, 

a calculer (2) pour chaque trajet les attenuations a„ 
et a calculer (3) la distance d separant la station 
emettrice/receptrice du radio-mobile en conside- 
rant seulement les retards et les attenuations du 
trajet direct et d'un multitrajet quelconque. 

Applications : Localisation et suivi de radio-mobiles. 
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Description 



[0001] La presente invention concerne un procede de localisation de radio-mobiles en milieu urbain. 
[0002] Elle s'applique notamment a la localisation et au suivi de radio-mobiles dans un reseau GSM pour anticiper 
s leur passage dune cellule a la suivante et ameliorer la reception des stations de base par la connaissance de la posi- 
tion des radio-mobiles emetteurs. 

[0003] Dans des configurations simples de propagation, la localisation de radio-mobiles est traditionnellement obte- 
nue a I'aide de radio-goniometres disposes en des lieux determines d'un territoire. La position d'un radio-mobile peut 
alors etre determinee par triangulation en reportant sur une carte geographique les lieux ou sont situes les radiogonio- 

io metres ainsi que les directions d'arrivee sur ceux ci des ondes radioelectriques emises par le radio-mobile. 

[0004] Cependant, cette methode qui repose sur Texploitation des seuls chemins de propagation directs apparatt dif- 
f icilement applicable dans le cas ou les ondes empruntent plusieurs trajets et est egaJement lourde de mise en oeuvre 
car elle necessite plusieurs stations de reception. Ce type de propagation est notamment celui qui est rencontre en 
milieu urbain dans les reseaux de communication radio-mobiles cellulaires, ou le spectre des communications est situe 

is dans une gamme d'ultra hautes frequences s'etendant typiquement de 30 MHz a 3 GHz pour le reseau GSM. Pour 
cette gamme de frequences la propagation des ondes radioelectriques entre un emetteur et un recepteur peut 
empreinter des trajets multiples dus aux reflexions et aux diffractions sur des obstacles interposes sur le chemin de pro- 
pagation direct entre I'emetteur et le recepteur. La gamme des frequences utilisees est aussi la cause de trajets diffus, 
aussi bien temporellement que spatialement, qui degradent les performances des liaisons. 

20 [0005] Le but de I'invention est de pallier les inconvenients precites. 

[0006] A cet effet, I'invention a pour objet un procede de localisation spatio-temporelle d'un radio-mobile relatrvement 
a une station de base emettrice/receptrice caracterise en ce qu'il consiste : 

a calculer la reponse impulsionnelle du canal de propagation entre le radio-mobile et la station emettrice/recep- 
25 trice, 

a estimer a partir de la reponse impulsionnelle les couples (angles, retard) du trajet de propagation direct et des 
multitrajets resultats de reflexions sur des obstacles situes dans le champ de propagation des ondes entre la sta- 
tion de base et le radio-mobile. 

a calculer I'attenuation a t de I'onde propagee par chaque trajet, 
30 - et a calculer la distance d separant la station emettrice/receptrice du radio-mobile en considerant seulement les 
retards et les attenuations du trajet direct et d'un trajet non direct quelconque. 

[0007] D'autres caracteristiques et avantages de I'invention apparaTtront a I'aide de la description qui suit faite en 
regard des dessins annexes qui represented : 

35 

la figure 1 le procede selon I'invention sous la forme d'un organigramme, 
la figure 2, un module de reception monocapteur, 

la figure 3, une illustration d'une reception multicapteur par une antenne circulaire, 

les figures 4 et 5, une illustration d une transmission a deux ou trois trajets en milieu urbain. 

40 

[0008] Le procede selon I'invention se deroule de la fagon representee a la figure 1 suivant trois etapes principales. 
Dans une premiere etape, representee en 1 a la figure 1 , il consiste a calculer la reponse impulsionnelle h(t) du canal 
de propagation pour estimer a I'etape 2 les angles 6 m et les retards - m des differents trajets ainsi que les attenuations 
a n des ondes propagees sur chacun des trajets et enf in a calculer a I'etape 3, a partir des retards x m calcules et de la 
45 reponse impulsionnelle h(t) du canal, la distance d separant un radio-mobile de la station d'emission/reception. 

[0009] Le calcul de la reponse impulsionnelle du canal qui est effectue a I'etape 1 a lieu sur le signal obtenu a la sortie 
d'une antenne multi-capteurs. 

[001 0] En presence de plusieurs trajets M d'emission bande etroite de radio-mobiles, le signal obtenu a la sortie d'un 
capteur n de I'antenne multi-capteurs est egal a la somme des M signaux s m (t) emis dans les directions e m retardes de 
so x m et multiplies par la reponse a t (Q m ) du capteur i dans les directions e m a laquelle s'ajoute un bruit bj(t). Dans I'hypo- 
these d'une transmission numerique a bande etroite de frequence, une representation de la sortie vectohelle de 
I'antenne composee de N capteurs s'obtient en concatenarrt la sortie xj(t) de chaque capteur, chaque sortie X|(t) veri- 
f iant la relation : 
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£s m (t).a i (6 m )+b i (t) 
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[001 1 ] En ihtroduisant un vecteur directionnel a(9 m ) tel que 

a<e m ) = [a,(e m ) a N (ej]' (2) 

et en concatenant la sortie Xj(t) de chaque capteur, on obtient un vecteur x(t) tel que : 

M 

x(t) = [x^t), Xi (t). x N (t)] 1 = Y. a(e m )s m (t) + b(t) (3) 
[001 2] En rassemblant les vecteurs directionnels a(6 m ) dans une matrice 

A = [a(e 1 ),a(9 2 ) a(9 N )] 

la relation precedente s'ecrit encore : 

x(t) = As(t)+b(t) (4) 



15 
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25 avecs(t) = [s 1 (t).s 2 (t). • -s M (t)]\ 
ou $1 est le signal du 1 er trajet. 
s 2 est le signal du 2'® me trajet, 



30 



s M est le signal du M idme trajet. 

[0013] Suivant le modele de propagation qui est represents schematiquernent k la figure 2, le signal x(t) qui est 
obtenu a la sortie d'un capteur de reception est un signal de la forme : 

M 

35 x(t) = £ a m .e'* ra .a(e m ).s(t - tJ + b(t) (5) 

ou designe I'attenuation du m* me trajet 
ao T m est le retard de groupe du m*"™* trajet 
M est le nombre de trajets 
et <j> m est le d6phasage du m ,6mo trajet. 



45 



soit encore sous la forme d'une convolution du signal emis s(t) avec la reponse impulsionnelle h(t) du canal : 

soitx n (t) = h n (t)*s(t) + b n (t) (6) 



ou b n (t) repr^sente I'enveloppe complexe du bruit adcfitrf. 
so [001 4] Dans le cas d'une reception sur I'antenne composee de multi-capteurs. comme cela est represents a la figure 
3. la reponse impulsionnelle du canal est modifiee par la reponse de I'antenne suivant les directions d'arrivSe. Cette 
reponse est de la forme : 

m 

55 h nW= £a n (e ro . Ym )a m e'* m 6(t-x m ) (7) 

m»1 
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en supposarrt le canal ideal, c'est-a-dire de bande irrf inie. 
[0015] En mode echantillonne. la reponse du canal de bande B prend la forme : 



h„M = F e £ a n (0 m )a m e' Vm sin c((k -F e xJ B) 



(8) 



avec k € Z, I'espace des nombres errtiers. 



et sin c((k- F e T m ) B) : 



Sin(rc(k-FeT m )B) 
"(k-F e x m )B 



(9) 



15 F e est la frequence d'echantillonnage 

B est la bande en frequence reduite du signal. 

[001 6] Cette reponse peut etre estimee par un f iltre causal, a reponse impulsionnelle f inie de longueur L et de coef- 
ficients tels que : 

20 



h„M =p n (9 m )a m e $m sin C ((k - F e x m ) B) 



(10) 
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35 



pour k € {0 ... L} 

et h n [k] = 0 pour k (0 ou k >L 



Le vecteur h n correspondant s'ecrit h n = [h n (0), • • • , h n (L)] x 



(11) 



[0017] Suivant un premier mode de realisation de I'invention, I'estimation des angles 6 m et des retards x m de Tetape 
2, est executee de fa?on disjointe, en utilisant la relation matricielle suivante : 



H = [h 1 .h 2 ,...,h N ] = IVI.aA t (12) 



avec : 
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M = 



sinc(-T-|F e B) sinc(-T 2 F e B) 
sinc((l-t 1 F e )B) sinc((l-x 2 F e )B) 



sinc(-T M F e B) 
sinc((l-T M F e )B) 



sinc^L-x^eJB) sinc((L - T 2 F e ) B) sinc((L- t M F e ) B)_ 



(13) 
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n = 



70 



ai 0 -•• 0 
0 a 2 : 
: 0 
0 ••• 0 ajsj 
A = [a 1 (9 1 )...a N (G N )] 



(14) 



(15) 



15 



[0018] L'estimation des parametres temporels x m s'obtient a partir de la relation (12) en estimant un projecteur 

P HHt 



sur lespace temporel HH t et en appliquant soit la methode connue sous Pabreviation anglo-saxonne MUSIC, 
de < <Multi Signal Classification >>en cherchant les maxima de la quantite 



20 



1 



(16) 
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m(T) = [sinc^-xFeB), .... sinc((L - xF e )B)]* 



35 



est une colonne de la matrice M(13). soit la methode deterministe connue sous le nom de < < maximum de vraisem- 
blance )) en utilisant un critere de la forme 

J(*1 T M )=t r ((l L+1 -M(MtM)- 1 M+)HH^) (17) 



40 [0019] Une description de la methode MUSIC peut §tre trouvee dans Particle de M. Ralph O. Schmitt ayant pour titre 
< ( Multiple source DF signal processing : An experimental system ) > publiee dans IEEE transaction on antennas and 
propagation AP34 N° 3 March 1986. 

[0020] La determination des M couples (6 m , x m ) permet de localiser les radio-mobiles. Dans le cas ou, comme cela 
est represents a la figure 4, seulement deux trajets, un trajet direct et un trajet comportant une reflexion sont selection- 
's nes, la determination des couples a lieu en considerant quatre parametres 0 1f , f^, Qui se ramenent a la connais- 
sance de 0^ , 6 2 et des retards t 1 . x 2 . En appelant d la distance du trajet direct qui separe I'gmetteur radio-mobile 4 d'un 
recepteur 5 et d 1t d 2 les distances du trajet indirect separant respectivemerrt le recepteur 5 d'un obstacle 6 sur lequel 
Ponde indirecte est ref lechie et Pemetteur 4 de I'obstacle 6. la distance d est fournie par la relation : 



so 



d= , 2 1 (18) 



K 1 

V a 2 V 



ou a<[ et a 2 designent les attenuations du premier et du deuxieme trajet et n la perte due k la reflexion du second 
trajet Les attenuations a n et ag sont obtenues a partir des couples (9 1t tj ; 6 2 . x 2 ) et de la relation matricielle 
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n=M # HA #t = 



"ai 0 1 
_0 a 2 J 



(19) 



5 



ou correspond a I'attenuation du i* me trajet 
[0021] Comme la perte de puissance dans une propagation libre se fait en 



10 



1 



15 



ou r est la distance parcourue, et c*2 sont lies par la relation : 



20 



*2 " d 2 



(20) 



avecd 



(T 2 -T 1 ) 2 .c 2 4>2d(T 2 -x 1 )c 
2[( T 2 ' T 1 )c+d-dcos( 6 2 -6 1 )] 



(21) 



25 



etd 2 = A /d 2 +d 2 - 2dd ^05(62 - 



(22) 



[0022] S'il y a plus de deux trajets, la figure 5 montre que Information apportee par tout trajet supplemental se 
30 resume a deux parametres, et !es nouvelles inconnues sont egalement au nombre de deux. En consequence, une 
detection de plus de deux trajets n'a pas d'interet dans le sens ou elle ne supprime pas les inconnues. II suffit dans ces 
conditions de limiter les traitements en considerant seulement le trajet direct et un trajet autre quelconque. Connaissant 
les angles d'arrivee des differents trajets, il est possible, gr§ce a une carte topographique detaillee en trois dimensions 
de trouver les points de reflexion les plus probables et ainsi de remonter a la source d'emission. 
35 [0023] D'autres modes d'estimation des parametres sont Egalement possibles. 

[0024] L'estimation disjointe des parametres decrrte prec6demment a I'avantage de reduire I'espace de recherche. 
Cependant, dans cette methode , les angles et les retards ne sont pas associes si bien qu'elle souffre d'un manque de 
resolution pour estimer les retards proches. 

[0025] D'autres modes de mise en oeuvre de I'invention permettent d'associer les parametres en estimant la matrice 
40 des attenuations il avec les differentes combinaisons (6 m , T m ) possibles en cherchant les couples rendant la matrice 
ft la plus diagonale possible. 

[0026] Selon un deuxifcme mode de mise en oeuvre de I'invention, le proc6de op6re par estimations successives en 
estimant successivement les angles 6 et les retards t tout en les associant. Pour cela, deux methodes peuvent §tre uti- 
lisees. La premiere consiste k estimer les parametres temporels x m connaissant les angles 6 m . 
45 [0027] En estimant des angles 6 m par la methode MUSIC ou la methode du maximum de vraisemblance l'estimation 
des retards x m est obtenue en estimant la matrice M par la relation matricielle 



50 




(23) 



[0028] Chaque colonne de la matrice M forme un projecteur permettant de trouver le retard correspondant en 
appliquant l algorithme MUSIC pour une seule source demission. 
55 [0029] Le critere utilise est de la forme : 
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P(t) = 



m^(t) rn(T) 
ID(ti)m t (ii) 



(26) 



70 



75 



D3 T (tj)m(T j ) 
avec m(x) tel que M = (m(Tj) ... m(-t M )] 

sinc(-xF e B) 
et m{x) = : 

sinc((L-TF e )B) 



20 



[0030] L'estimation des retards peut egalement etre obtenue par correlation. En considerarrt que le signal recu sur 
un capteur de M sources est de ia forme : 



x(t)o £a(e m )s m (t) + b(t) 

m-1 



(25) 



25 



30 



et connaissant les angles 6 m . il est possible de considerer un signal : 



s[k]= A # x[k] 



qui est forme de M versions du signal emis retardees tel que : 



35 



AO 



s(k-t 2 F e ) 
s(k-T M F e ) 



(26) 



[0031] ^autocorrelation d'une des composantes de s [k] permet d'obtenir une reference temporelle et I'intercorrela- 
tion de cette composante avec les autres composantes donne une estimation des retards relatifs. 
45 [0032] De facon similaire a la methode precedente. il est aussi possible d'estimer les parametres spatiaux 0 m con- 
naissant les retards t m . Ceux-ci sont obtenus en calculant Pequation matricielle : 



so 



(27) 



[0033] Chacune des methodes utiiisees dans ce deuxieme mode de realisation du precede selon Pinvention permet 
apres avoir estime les angles ou les retards de trouver les autres inconnues par M recherches en une dimension. Elles 
ont pour avantage de repousser certaines ambiguTtes, mais la seconde recherche apparalt tres dependante de la qua- 
lite de la premiere estimation. Une aHernative serait d'effectuer, comme cela a ete decrit par M. J. GROUFFAUD dans 
son memoire de these de Doctorat de TEcole Normale Superieure de Cachan de juin 1997 intitulee <( Identification 
spatio-temporelle de canaux de propagation a trajets multiples > ) . des recherches successives en alternant les metho- 



7 



10 



15 



20 



30 



EP0 926 510 A1 

des pour affinerla precision de I'estimation des parametres. 

[0034] Enfin, suivant un troisieme mode de realisation du procede selon Invention, I'estimation des parametres e m 
et x m a lieu de tegon conjointe en appliquant des methodes de determination des reponses impulsionnelles fondees sur 
une decomposition en sous-espaces, un sous-espace signal de dimension M egale le nombre de sources et un sous- 
espace bruit, complementaire du sous-espace signal de dimension N(L + 1)-M, ou en appliquant la methode du maxi- 
mum de vraisemblance. 

[0035] Pour estimer conjointement les angles et les retards par la methode des sous-espaces, un projecteur P Rh est 
cree a partir de la matrice d'autocorrelation des firtres 



[0036] (.'estimation b du vecteur h = [hJO], h,[1] h^L] h N [0]. h N [1] h N [L]]' a lieu en recherchant les 

maxima de la quantite 

POT = gt(9,) pLs(9,) <28) 



ou g(6.T) = a(e)«)m(T) et 

25 M 



m=1 



avec 



[0037] Dans la recherche des parametres par la methode du maximum de vraisemblance les parametres cherches 
35 sont les minimum du critere, 

ouG = [g(0 1 ,T 1 )...g(9 mi T m )] 



45 [0038] Cette recherche a lieu dans Tespace de dimension M x 2 des parametres (e 1t t 1 , • • • , 0^, x M ). 

[0039] L'optimisation des drfferents modes de realisation du procede selon invention peut avoir lieu en appliquant 
une methode connue d 'optimisation du type Gauss-Newton. Apres avoir trouve un point d'initialisation par les methodes 
decrites precedentes, une optimisation des angles e m peut avoir lieu en deroulant I'algorithme suivant : 

so Initialisation : e 0 = e, nit 

Iteration k : 

Calcul du gradient v k-1 et du hessien H k ., au point 6 k .. 
Determiner 

55 A ~ -1 
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Test de la convergence 

si O^ 7 k -iB ( e STOP sinon k = k + 1 

c 

[0040] Le meme type d'optimisation peut etre utilise pour optimiser la recherche des retards t m . 
[0041] Un autre mode d'optimisation peut etre mis en oeuvre par une version approchee du maximum de vraisem- 
blable suivant par exemple la methode decrite dans IEEE transaction on signal processing, Vol. 39, N1 1, Novembre 
1991 des auteurs M. Viberg. B. Ottersten, T. Kailath ayant pour titre : < ( Detection and estimation in sensor arrays using 
io weighted subspace fitting )>. 
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ANNEXE A 
GLOSSAIRE 

A.1. Conventions standards : 

x* : conjuguee de x 

x ! : transposee du vecteur x 

x* : transposee conjuguee du vecteur x 

X* : transposee de la matrice X 

X 1 * : transposee conjuguee de la matrice X 

X" 1 : inverse de X 

X* : pseudo-inverse de X : X # = (x t x)" 1 X t 

tr(X) : trace de X 

E[.] : esperance mathematique 

x, X : valeurs estimees 

x* y : convolution des signaux x et y 

5(t) : distribution de Dirac 

sinc() : fonction sinuscardinal 

® : produit de Kronecker 

A.2. Constantes entieres : 

N :nombre de capteurs composant I'antenne 

M : nombre de trajets 

L : longueur de la reponse impulsionnelle 

K : nombre d'echantillons temporels pour composer I'observation 
spatio-temporelle. 

A.3. Variables scalaires : 

x n (.) observation sur le n* me capteur 
s() : signal emis par la source 
b() : bruit 

a 2 : variance du bruit 

h n () : reponse impulsionnelle du n* mo canal 
F. : frequence d'echantillonnage 
T 0 : periode d'echantillonnage 
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a m : amplitude du m i6mtt trajet 

<J> m : dephasage du m ,6mo trajet 

6 m : azimut du m* me trajet 

y m : site du m i6me trajet 

x m : retard de groupe du m ,emc trajet 

a n (0) : gain du n* me capteur dans la direction (9) 

J() : critere a optimiser 

A.4. Variables vectorielles et matricielles : 

x :observation spatiale sur I'antenne, de longueur N 

y observation spatio-temporelle sur I'antenne, de longueur NK 

h : vecteurdes reponses impulsionnelles, de longueur N(L+1) 

a(9) : vecteur directionnel spatial, de longueur N 

m(x) : vecteur tempore!, de longueur L+1 

l q : matrice identite de taille qxq 

R x : matrice de covariance du vecteur x 

R y : matrice de covariance du vecteur y 

Rh : matrice de covariance du vecteur h 

M : matrice des M vecteurs temporels : M = [m(Ti), — , m(^M)] 

A : matrice des M vecteurs directionnels : A = [a{Q^) t — , §(9^)] 

n : matrice des attenuations des M trajets 

H : matrice des filtres mis en colonnes, de taille (L+1)xN 

G : matrice formee des reponses impulsionnelles 
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ANNEXE B 

LE MAXIMUM DE VRAISEMBLANCE DETERMINISTE 

5 

En considerant un signal de la forme : 
y k = A(0)s k + b k 

10 et le bruit comme etant gaussien, centre de variance a 2 , la 

distribution de probability du bruit est de la forme : 

FT 

La distribution de probability des donnees qui est a maximiser est 
donnee par la relation : 

20 „ I x 1 -^K-ASkI 2 

f(Xk|Sk) = ^ 6XP 

La vraisemblance des donnees pour le vecteur des donnees 
25 concatenees est alors : 

K 1 -^|x k -As k | 2 

L ( x 1- X 2 x k) = n-j^fT2j eXP ° 

k-ut M |o' e | 

30 

Pour maximiser la vraisemblance cela est equivalent a minimiser 
la relation : 

-logfLfxL x 2f .... x k )) = -MKIoga 2 + f> k - As k | 2 
35 k=1 

La minimisation de ce critere par rapport a 0 amene a minimiser : 

K 2 

*K = Zl x *- As k| 
40 k = l 
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d'ou : 



k=U 
comme 



|l-AjA f Aj"V 
|u k | 2 = uju k = tr(uju k ) = tr(u k u - ) 



ZJ U M 
k=i 
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eit devient : 
K 



e k = Z tr ( x M n bX k ) = tr(n b R x ) 



Revendications 

1. Procede de localisation spatio-temporelle d'un radio-mobile relativement a une station de base emettrice/recep- 
75 trice caracterise en ce qu'il consiste : 

a calculer (1) la reponse impulsionnelle du canal de propagation entre le radio-mobile et la station emet- 
trice/receptrice. 

a estimer (2) a partir de la reponse impulsionnelle les couples (angles, retard) du trajet de propagation direct 
et des multitrajets resultats de reflexions sur des obstacles situes dans le champ de propagation des ondes 
entre la station de base et le radio-mobile, 
a calculer (2) pour chaque trajet les attenuations a,, 

et a calculer (3) la distance d separant la station emettrice/receptrice du radio-mobile en considerant seule- 
ment les retards et les attenuations du trajet direct et d'un multitrajet quelconque. 

Procede selon la revendication 1 caracterise en ce qu'il consiste a determiner la distance d par la relation : 

d _ (^2- T i)c 

ou x 2 - ^ est I'ecart des temps de propagation entre le trajet direct et le trajet refiechi, 
c est la celerite des ondes transmises sur le canal de propagation, 

35 

ill 
a 2 



20 



25 



AO 

est le rapport des attenuations des propagations sur le trajet direct et le second trajet, 
et u la partie due a la reflexion de I'onde sur le second trajet. 

3. Procede selon la revendication 2 caracterise en ce qu'il consiste pour calculer (2) les attenuations et a 2 des 
45 ondes sur le trajet a calculer une matrice diagonale Q def inie par le produit matriciel : 

50 OU 

55 

est la matrice des coefficients de la reponse impulsionnelle de canal, 
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H = [h 




et A 1 est la matrice transposee des vecteurs directionnels (e n , y n ) correspondant aux multitrajets. 



4. Procede selon Tune quelconque des revendications 1 a 3, caracterise en ce qu'il consiste a determiner ia reponse 
impulsionnelle du canal (1) au moyen dune antenne multi-capteurs coupe a un filtre causal, a reponse impulsion- 



5. Procede selon Tune quelconque des revendications 1 a 4, caracterise en ce qu'il consiste a effectuer une estima- 
tion des angles et des retards en estimant les parametres temporels a partir de la connaissance des angles e m . 

6. Procede selon I'une quelconque des revendications 1 a 4, caracterise en ce qull consiste a effectuer une estima- 
tion des angles et des retards en estimant les parametres spatiaux 6 m a partir des retards i m . 

7. Procede selon Tune quelconque des revendications 1 a 4, caracterise en ce qu'il consiste a effectuer une estima- 
tion conjointe des angles 0 m et des retards T m . 

8. Procede selon I'une quelconque des revendications 1 a 7, caracterise en ce qu'il consiste a estimer les 
parametres temporels et les directions d'arrivee des multitrajets par la mise en oeuvre de I'algorithme connu sous 
le nomde << MUSIC >). 

9. Procede selon Tune quelconque des revendications 1 a 7. caracterise en ce qull consiste a estimer les parametres 
temporels et les directions d'arrivee des multitrajets par la mise en oeuvre de I'algorithme deterministe connu sous 
la denomination de < < maximum de vraisemblance ) >. 



nelle finie. 
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